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Abstract 

The present study examines the impact of large- 
scale cross flow on the creation of ice roughness ele- 
ments on the leading edge of a swept wing under 
glaze icing conditions. A three-dimensional triple- 
deck structure is developed to describe the local in- 
teraction of a 3D air boundary layer with ice sheets 
and liquid films. A linear stability analysis is pre- 
sented here. It is found that, as the sweep angle in- 
creases, the local icing instabilities enhance and the 
most linearly unstable modes are strictly three di- 
mensional. 

Nomenclature 

T m bulk freezing temperature, 273°K for water 

y si solid-liquid surface tension 

Y iv liquid-vapor surface tension 

H 5 / latent heat of freezing per unit mass of ice 

H/ v latent heat of evaporation per unit mass of water 

q density 

c p specific heat 

k thermal conductivity 

a thermal diffusivity 

p dynamic viscosity 

0 m density ratio, = 1 - (qi/q s ) 

LWC liquid water content of impinging droplets 

|3(s) local collection efficiency 

x sum of the principal curvatures of the interface 

Uqo ambient airspeed 

Tqo ambient air temperature and airspeed 

Poo ambient air pressure 

L re f characteristic length scale 

AT freestream supercooling, = T m - Too 

Mqo freestream Mach number 

Pr Prandtl number of air, = 0.72 

Re freestream Reynolds number, = QooUooL/goo 

Ste Stefan number, - C pyS AT/H S/ 

Sc Schmidt number, = \i/qDab 

Dab mass diffusivity coefficient 

M relative molecular mass 

JWd viscosity ratio of water to air 

P v saturated vapor pressure of pure water 


(U,T) e inviscid surface velocity / temperature 
(s,N) boundary layer scaled coordinates 
(u,v)b boundary layer scaled velocity components 
(P,T)b boundary layer scaled pressure / temperature 
ho(s,t) boundary layer scaled water film thickness 
6o(s,t) boundary layer scaled ice sheet thickness 
(X,Y) local triple deck scaled coordinates 
(U,V) triple deck scaled velocity components 
(P,T) triple deck scaled pressure / temperature 
h(X,t) triple deck scaled water film thickness 
6(X,t) triple deck scaled ice sheet thickness 
Vimp scaled impingement velocity 
5 W all boundary layer scaled wall layer thickness 
T w ,s scaled steady state inner wall temperature 
to dimensionless global time 
t dimensionless local triple deck time 
d roughness element diameter 
R airfoil leading-edge nose radius 
A c accumulation parameter, dimensionless 
n freezing fraction, dimensionless 

Subscripts 

oo in the freestream air 

/ in the liquid water region 

s in the frozen ice layer 
w in the airfoil surface metal wall region 

Superscript 

* dimensional quantity 

Introduction 

This work analyzes the linear stability of a three-di- 
mensional laminar boundary layer flowing past a 
thin layer of surface water, formed and replenished 
from the impingement of small supercooled drop- 
lets, on a growing glaze ice substrate near the lead- 
ing edge of a wing surface with the presence of sig- 
nificant cross flow. Specifically flow over a swept 
wing is considered, as shown in Fig. 1. Recent ex- 
perimental investigations on glaze ice shape forma- 
tion on swept wings 1_4 have shown that the sweep 
of a wing can affect the glaze ice accretion process. 
Figures 2 and 3 are front-view images of a portion 
of the wing leading edge after approximately 1 min- 
ute of ice accretion under similar icing conditions, 
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with n=0.41 and 0.43 respectively. Clearly the in- 
duced cross flow on swept wings has changed and 
modified the surface water run-back behavior, rela- 
tive to the zero sweep configuration, and conse- 
quently the glaze ice accretion process is changed 
as well. One of the most noticable difference be- 
tween Figs. 2 and 3 is the replacement of a smooth 
zone (i.e. 2D dominant modes with diminishing 
growth rate) with well defined 3D roughness field 
along the attachment line (or the stagnation line) re- 
gion when the wing is swept. 

In a previous study on the formation of ice rough- 
ness on a NACA 0012 airfoil leading edge under 
glaze icing conditions, Tsao & Rothmayer 5 showed 
that a mutual interaction between the film dynamics 
and the ice accretion process can lead to the creation 
of surface ice roughness elements. The linear stabil- 
ity analysis has predicted that the most unstable and 
fastest growing modes have streamwise wave- 
length of millimeter size, and also there is a narrow 
region of strong instability centered at some dis- 
tance downstream of the nose. Although qualitative 
features of glaze ice roughness formation on an air- 
foil leading edge observed experimentally 6-7 are 
well captured by the 2D model, its usefulness in 
helping one to understand the aforementioned 
change in roughness formation mechanism associ- 
ated with swept wings is then limited because the 
flow is 3D, even if only weakly so. Therefore, a ful- 
ly three-dimensional triple-deck model is used to 
study the effect of cross flow on local icing instabi- 
lities leading to the formation of glaze ice rough- 
ness. 

The second area of interest of the paper concerns the 
possible connection between the Messinger freez- 
ing fraction and the icing instability. The Messinger 
model in its original form represents an equilibrium 
energy balance, in which the freezing fraction indi- 
cates how much super-cooled water droplets upon 
impact will turn into ice mass in the equilibrium 
steady state. Its numerical value, ranging from 0 to 
1, determine the physical appearance of the accreted 
ice. But why does this parameter have such a mor- 
phological feature? Is it related to the dynamic as- 
pect of the icing process? As a starting effort, a 
theoretical derivation of the freezing fraction ex- 
pression within the context of the current model is 
made , and it is found that the freezing fraction is di- 
rectly tied to the samll scale icing instability. A brief 
discussion of its implication will be given. 
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The Nonlinear 3D Triple-Deck Flow 

Tsao & Rothmayer in Ref. 5 have laid out a com- 
plete derivation of the overall asymptotic structures 
for a 2D airfoil glaze icing problem, which includes 
global boundar-layer and local triple-deck scaled 
structures. As for the 3D counterpart problem, i.e. 
flow over a swept wing, the asymptotic scaling 
analysis given in Ref. 5 can be readily applied. Un- 
fortunately, space does not permit a full discussion 
and so the paper will focus only on the local flow 
structures which are important to the glaze ice 
roughness formation process. 

In describing the local flow, a Cartesian coordinate 
system (x,y,z) and their corresponding velocities 
(u,v,w) are taken to be centered at some point on the 
wing surface, see Fig. 1. The local near-wall 
streamline (for y— ►()) defines the streamwise direc- 
tion, x, to which the crossflow direction, z is nor- 
mal, and y denotes the coordinate which is perpen- 
dicular to and out of the wing surface. Notice that 
the angle between the local streamwise direction 
and the global chordwise direction is the flow angle 
0, where 

0 = tan -1 (U s /U c ). 

Here the subscripts c and s indicate the chordwise 
and spanwise directions respectively. Also, it 
should be noted that the present analysis neglects 
possible surface curvature effects, which is general- 
ly acceptable in the thin layer model considered 
here but may fail in other situations. 

As shown in Fig. 4, the local air flow is a steady 
three-dimensional tripl-deck structure. The basic 
equations and boundary conditions governing the 
local 3D triple-deck flow is similar to the 2D case 
derived in Ref. 5 (i.e. the problem in the z-direction 
is essentially a mirror image of the problem in the 
x-direction provided finite boundary-layer scaled 
wall shear stressess and edge velocities) and will 
simply be stated here. Interested readers are referred 
to that study, and some related detailsof 3D bound- 
ary layer instability problem may be found in Stew- 
art & Smith 8 and Ryzhov & Terent’ev 9 . 

Let us further confine ourselves to the lower near- 
wall viscous sublayer ( lower-deck ) where the air 
flow can affect the growth dynamics of the water 



film and ice layer. The governing equations for the 
local triple-deck structure are given in Prandtl 
transposed variables by 

U x + V Y + W z = 0, (1) 

( 2 ) 

P 7 + Wyyj (3) 


UU X + VU Y + WU Z = - P x + Uyy, 


UW X + VW Y + WW Z = - 
and 

UT X + VT y + WT Z - 1 1 ± yy 

with the boundary conditions 
U(X,0,Z) = V(X,0,Z) = W(X,0,Z) - 0, (5) 


z 

Pr _1 T 


( 4 ) 


T(X, 0, Z) - 


+ Hv ; „ 


X 


^(T Y (X,0,Z) - E) . ... |mp 

(h(X, Z, t) - 6(X, Z, t)) + r(6 xx + 6 ZZ ) (6) 
and the farfield matching conditions 
U -»• X x (Y + h — h 0 + A(X, Z)) 

W \ z (Y + h - h 0 + A(X, Z)) 

T^q[Y + h-h 0 + A(X, Z) 


+ T^© (h 0 - 6o)], 


as Y-> oo, 
as Y -» oo, 

(7,8,9) 

as Y oo 


where f.\ and /./ denote the local boundary layer 
wall shears, q is the local boundary layer surface 
heat transfer, ( ho, 6o) represent the boundary layer 
scaled water film and ice sheet thickness and 


0 


1 + 


k Z H 0 v im P ,0 


kooq 


E 


c - n „ - 1 /2 H/v w evap L ref 

E = Re k„ AT ’ 


^evap 


h c /PrWM 


T [ 


- v,s 


H r (2f)P v ,c 


Sc/ \ M a / L p; 


H 0 = Re “ 2//8 Jb “ 3 Pr (^w-)[y— ^ 

0 v Poo V L 2c p ,ooAT C P ,°° J 


(v imp ,Veva P )o = Re 6 />l4E (^% 


m 


evap 


Pz '’p z Uc 


X 


H = Re 1 / 8 ^- 1 Pr(^^)[^T- - 

v Poo k/ L 2c p ,ooAT C P ,°° J 


(Vj 


m 


imp’ V evap, 


evap 




)• 


The viscosity ratio of water to air, Jfb, is assumed 
to be large so that the lubrication approximation can 
be used to describe the film flow. 


The unknown displacement function, A(X,Z), is re- 
lated to the self-induced air pressure through the 
Cauchy-Hilbert integral: 


P(X,Z) 


* 

_ - 1 
2jt 

- 


(U e a/a^ + W e d/dri) 2 d^dr) 
[(X - If + (Z - T]) 2 ] 1 / 2 


. ( 10 ) 


The U e and W e stand for the velocities just outside 
the boundary layer. The Gibbs-Thomson parameter 
r, shown in eqn. (22), comes from the Gibbs- 
Thomson relation and is defined as 

r r r>^1/4 a 

Q 0 

Ps H sZ AT L ref L ref 
where do is a capillary length defined by 
YszTmRe 1 /* 


= A (11) 


d 0 = 


Ps H sZ AT 


( 12 ) 


The local water film motion and ice layer growth 
are nonlinearly coupled and governed by 


77 , d 

ht + ax 


+ -77 
az 


~ 2 ~ 3 

U Y (X,0,Z)^- (P x + G x sin0 x )|- 


~ 2 ~ 3 

W Y (X,0,Z)|- - (P z + G z sin0 z )^ 


.(13) 


and 


Mmp v evap 


Ctoc 

ou 


\«56 

= 

ri -T w , s \ 

/ at 

k s \ 

( ‘Vail / 


k z/k 


^I^(T y (X,0,Z) - E) 


+ H V; 


imp 


(14) 


where h = h — 0 m 6 , h = h — 6 and 
P x = Px — T(Ah) x + G y (cos0 x + cos 0 z )h x 
P z = P z — T(Ah) z + G y (cos0 x + cos 0 z )h z 


with 


Ah — h X x + ^zz* 


The gravitational and surface tension parameters of 
liquid water are defined as 

(G X;Z ,G y ) = (Re-V 8 ,R e - 3 / 8 ) ( p zX ^) 

= (Re -17 / 8 , Re _19 / 8 )(p z poogL 3 /p^) 


and 


T = Re- |3 / 8 ( PooYzv L/pl ). 


Also in the Stefan equation, i.e. eqn. (14), the tri- 
ple-deck scaled Stefan number is defined as 
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where 


S 1 = Re^pHoSte)- 1 , 

_ C p,.s AT _ P 00 U 00 L ref 

- Ke “ ^ 


A = A [q0 + 



k 5 aooSjPr 
kooOs 0 


A = (1 -^)Ai'(O) 

Jy oo 


<P 2 

=A- Ai(0). 
0k 


Linear Stability Analysis 

As in the studies of Smith 10 and others 5 ’ 11 ’ 12 , the 
flow is resolved into its primary and perturbation 
parts. Although we are considering only the leading 
order term, Smith 10 shows that high order terms 
could be added to improve the accuracy of the solu- 
tion. The variables are expanded in normal mode 
form as follows: 

(U, V, W, P, P, T, A, h, 5) = (X x Y,0,X z Y,0,0, 
q[Y + ^0(h o - 6o)],0,h 0 ,6 0 + e(U(Y), (15) 

V(Y), W(Y), P, %), T(Y), A, h, &)e ! (k ' x + k ' Z) + wt , 

where k x and k z (both are real) denote the triple- 
deck scaled wave numbers in X and Z directions and 
(/*)x,z = 2 jt/(£ x ,£ z )* = 2 jtLtd/(^x^z) are the associ- 
ated disturbance wavelengths (here Ltd = 
Re -3/8 L re f). The triple-deck scaled frequency co is 
complex, co = co r + i co;, and may be readily related 
to the well-known complex wave speed, c = c r + i 
q , where c = co /(-ik) with k 2 = £ x 2 + k z 2 . As usual, 
c r is the phase speed of the wave, whereas q deter- 
mines the degree of damping (if q < 0, co r < 0), or 
amplification (if q > 0, co r > 0). 

The three-dimensional linear stability problem can 
be further reduced to an equivalent two-dimension- 
al problem if one lets 

k 2 = k K 2 + k z 2 , kU = k x U + k z W, 

k\ — & x X x + k^kjj k\Jt = k x U e + £ z W e 

cos0 = cos0 x + cos0 z , (16) 

£Gsin0 = £ x G x sin0 x + k z G z sin0 z , 

V = V, T = T, and P = P. 

The ensuing quadratic dispersion relation for the 
linear instability eigenvalues, i.e. co, is found to be 

(Y P l o) 0)2 + ( A P - 0mN o L o + e,yL 0 ) (0 

+ (Ae, - N 0 L 0 e 2 - M v ) = 0, (17) 


where 


n _ e 1Jt / 6 k 1 / 3 \k\ Ue 2 
1 3 X 5 / 3 Ai'(0) 


L 0 = L x Ai(0) - L 2 Ai'(O), 
k, <E> 9 

L 1 = TT“ + T - ( h 0 ~~ 6 o)> 

JY 00 JY 00 

L 2 = (l-p-)^(h 0 -b 0 ). 


Nn 




Oo 


Yt'i Yfc'o 

(1 - ^)Nj - =-^N 


0k c 

% 

r 


2 ’ 


N 1 = ai Ai'(0) + Pr4Ac,Ai(0) + PrlaL-(O) 


+ (?) (PL -T 11 c l L k (0), 

k p r3 


N 0 = 


= a,Ai(0) + Pr3aLj(0)f 0 


q,(Pr - 1), 


Prf 


c iLk(0)> 


Oi 


'Cn.coAT^Scj ImT 


3 /M w \P v >s Hr(_.«)Pv.. 


I 


d> 2 = 


q H /vVPrV/M w \ r 17.15 P v ,s AT 




c p ,ooAT \Sc) \M a ) L p; (T m - 38.25) 


d 


a! = L" 1 ! q0(3 + + Pr5L-(0) ct + 

k Prl 

q( pr - l)L k (0) _i] _ Li[pr i L(0) _ + 
k Prj 

q(Pr - l)L k (0) 

k Prf Cl]> ’ 

= e 3 (3 + ^i£(h 0 - 6 0 ) 2 0 A;(0) c, + 
±£ 2 (h 0 -6 0 ) 3 P, 


e 2 = i&k(h 0 — 6 0 ) — i&(h 0 — 6 0 ) 2 G sin 0, 
e 3 = i kk (h 0 — 6 0 ) + ^& 2 (h 0 — 6 0 ) 3 G Y cos0 
- i£(h 0 - 6 0 ) 2 G sin 0 + ±£ 4 (h 0 - 6 0 ) 3 T, 

M v = N 0 [AjAi'(0) + B m Ai(O)][q0 + ^Tk 2 ], 
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_ - J\p / P°° s./Pr\ 3 /M w \ p y.s H r ( p ,^)P vce 

PrR e V8Vp ; ^Scj \M a ) 1 P e * 


R ~ Jb r P-JPr\^M w \ 17.15 P v , s AT 
m Pr Re 1 / 8 ^ P/ A Sc/ \ M J l P; (T m - 38.25) J ’ 

with 


_ e 1IT / 6 A: 1 / 3 |A:| Ue 2 * . . ~ 

1 X2/3Ai'(0) ’ P= W Ue ' 

a = (qA) «» /6 *‘gWUe 2 


and 


0 = (y&^/V*/ 6 , 0 = (kk Pr) 1 / 3 e i31 / 6 . 


It is noted that, in the above expressions, Ai denotes 
the Airy function of the first kind. Two special func- 
tions, Lj and L k , are introduced and they satisfy 

Lj — i Lj = - 1, L/oo) = Lj(oo) = 0, 
and 

Lk - ?L k = Ai (Pr-^), L k (oo) = L k (oo) = 0. 

Their numerical solutions at ^=0 are given as fol- 
lows: 

Lj(0) = 0.643949, L](0) = 0.469446, 

L k (0) = 0, L k (0) = 0.169538, 

along with 

Ai(0) = 0.35502, Ai'(0) = - 0.25881. 

A number of parametric studies, based on the test 
condition given for Fig. 2, are performed to assess 
the cross flow effects, i.e. via varying W e and Xz, on 
instability. Several important findings are obtained 
from the results: 

(1) There are two types of instabilities, wave-like 
(cjl>i) and broad-band ( 0 ) 2 ) modes, both present in 
the flowfield(see Fig. 5). 

( 2 ) The most linearly unstable modes are strictly 
three-dimensional when the cross flow is present 
(see Fig. 6). 

( 3 ) The instabilities are enhanced with increasing 
growth rates and wave speeds as the cross flow be- 
comes stronger, see Fig. 7. 

( 4 ) There is a considerable change in the most ob- 
servable instability growing patterns, in terms of 
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the characteristic length scales, when the cross flow 
is non-zero, see Figs. 8 and 9. 

Comparison with Experiments 

As mentioned early in the introduction, this study 
was partly motivated by the aim to understand the 
glaze ice roughness formation process on swept 
wings. To apply the stability analysis, i.e. eqn. (17), 
to the experimental data of Fig.3, a “ quasi 3D ” syn- 
thetic global flow profile was constructed for points 
along a given chordwise cross section (i.e. a NACA 
0012 airfoil) of a swept wing, see Fig. 10. The com- 
bined 3D flow mainly contains a 2D boundary layer 
solution for flow past the parabolic leading edge of 
a NACA 0012 airfoil at zero angle of attack with 
some prescribed spanwise velocity distributions. 

Two important observations are noticed from the st- 
ability analysis. They are summarized briefly as fol- 
lows: 

(1) The leading edge region is covered entirely with 
three-dimensional ice roughness elements of dif- 
ferent sizes and aspect ratios, see Fig. 11. The “av- 
erage” characteristic length of the 3D roughness is 
about 1 mm. The smooth zone observed in Ref. 7 
(see Fig. 2) for zero sweep angle condition has dis- 
appeared. 

(2) There is a narrow region of strong instability 
centered at some distance ( ^8mm ) downstream of 
the nose, see Fig. 12. The roughness elements in 
this zone are expected to grow much faster and larg- 
er than roughness elements in other locations. In 
comparison with our prediction, the critical dis- 
tance (measured from the attachment line to the be- 
ginning of larger glaze ice feather zone) given by 
Vargas et. aim Ref. 1 is about 6 mm after 1 min of 
ice accretion. 


Connection with the Freezing Fraction 


In icing scaling, the freezing fraction parameter n, 
derived from the Messinger energy-balance analy- 
sis for freezing at an unheated surface, can be writ- 
ten in the form 13-14 , 




(18) 


The terms involved are cj), the water droplet energy 
transfer parameter; 0 , the air energy transfer param- 
eter; and b, the relative heat factor which was first 



introduced by Tribus et. al . 15 The analytical forms 
of these terms are not presented here. Interested 
readers are referred to Refs. 13 andl4. For more re- 
cent information, readers are referred to Ref. 16 as 
well. 


Since Messinger’s energy-balance analysis does 
not consider heat loss from water runback or to con- 
duction into the cold wall, ice substrate and water 
film, the freezing fraction parameter is effectively 
related to the overall heat transfer parameter of air / 
droplet mixture, q0, of the current model. After 
some working, the freezing fraction can be rewrit- 
ten as 


n = (SoPr)- 1 


a z koo\/ c p ,/\/ - q©\ 

a oo k z y y c p,sy y v i mp ,o J 


(19) 


A numerical verification of eqn. (19) is carried out 
by comparing the calculated freezing fraction val- 
ues at the stagnation point for (laminar) flow over 
a NACA 0012 airfoil against values given by Refs. 
17 and 18 using eqn. (18). The results, given in Fig. 
13, show good agreement consistent over all test 
conditions. 


With this new form for q0, the large scale ice sub- 
strate growth rate at a given location on the surface 
can be rewritten as (see Ref. 5 for derivation) 


PrS 


o 


/aoo\ dSp 

\a s / at 




( 20 ) 

n 


where 


n = Reie JWj2 I = 



Pr S qCX oo 

/jia s L ref U c 


provided that 

P = (T m — T w (s)) 1? 

0 = (0w c p,wk w /0s kj) 1 / 2 P- 


The first term on the r.h.s. of eqn. (20) is the tran- 
sient growth from the cold wall, and the second 
term is the steady-state growth rate, i.e. as t — >oo. 
This suggests that the Messinger model really is a 
large time asymptote of the current model. Also the 
dispersion relation of eqn. (17) can be approxi- 
mated as 


C0j ~ => (X,Ue)~T 

^ «l e >" r 


via an order of magnitude analysis of the equation. 
With this expression for 002 , the freezing fraction is 
directly tied to the broad-band icing mode. Fur- 
thermore, by setting 002 = 0, the neutral stability 
wave number £2 is found to be 


^2 






Prv„ 


imp,0 


Re~'/ 2 L ref \ 

d 0 ) 


which gives 

h ~ ~ 0<nm) 


the smallest unstable ice roughness a physical 
length scale that is close to an ice feather. In summa- 
ry, the freezing fraction is found to be an important 
parameter that can control small scale ice roughness 
(feather) growth. 


Conclusions 

A three-dimensional triple-deck local flow struc- 
ture is developed to study the impact of large-scale 
cross flow on the creation of ice roughness elements 
on the leading edge of a swept wing under glaze ic- 
ing conditions. A linear stability analysis is per- 
formed locally and found that when the cross flow 
is present, the instabilities enhance and the most 
linearly unstable growing ice roughness are fully 
three dimensional. 

For glaze icing conditions tested on a NACA 0012 
swept wing 1-4 , two important observations are ob- 
tained from the stability analysis. First of all, the 
leading edge region is covered entirely with 3D ice 
roughness of mm size. The smooth zone, observed 
and reported by Anderson 7 , has disappeared. Sec- 
ondly, there is a narrow region of strong instability 
located at about 8 mm away from the nose. The r- 
oughness elements in this zone are expected to grow 
much faster and larger than roughness elements in 
other locations. In comparison with our prediction, 
the critical distance (measured from the attachment 
line to the beginning of larger glaze ice feather zone) 
given by Vargas et. al.m Ref. 1 is about 6 mm after 
1 min of ice accretion. 
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A connection between the Messinger freezing frac- 
tion and the current glaze icing model has been es- 
tablished. This parameter is found to be important 
for small scale ice roughness (ice feather) growth. 
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Figure 1: Schematic diagram of the flow geometry 
and the associated Cartesian coordinate system in 
which the y-direction is out of the plane. The bro- 
ken lines represent the surface streamlines. 



Figure 2: A close-up front view of the leading edge 
of a wing section, at zero sweep angle, after 1 min 
of ice accretion showing a smooth zone centered at 
the stagnation line and surrounding rough zones 
with roughness elements of mm size. 6-10-96 Run 
6 with n=0.41, taken from Ref. 7. 



Figure 3: A close-up view of the leading edge of a 
wing section, at 45° sweep angle, after 1 min of ice 
accretion showing roughness elements of mm size 
covering the entire leading edge region. 6-18-96 
Run 3 with n=0.43, taken from Ref. 1. 


Y, V 



Figure 4: Schematic diagram of the local 3D triple- 
deck structure with an oncoming steady3D bound- 
ary layer base flow. The (0x?6z) are the surface in- 
clination angles w.r.t. the horizon at a given point 
along the X and Z directions. 
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Figure 5: Typical wavelike (mode 1) and broad- 
band (mode 2) instability modes for a linearly un- 
stable case. 
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Figure 7:The most linearly unstable modes at dif- 
ferent water film thickness for (a) smaller cross flow 
(W e =0.26, kz=0.1)(b) larger cross flow (W e =0.7, 
kz=0.5) conditions. 
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Figure 6: The most linearly unstable modes at dif- Figure 8: The most linearly unstable wave-like 
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(a) 


(b) 


(c) 


Figure 9: The most linearly unstable broad-band- 
mode at different film thickness for (a)W e =0.26, 
X z =0.1(b) W e =0, X z =1.0(c)W e =5.0, X z =5.0 con- 
ditions. 


a) 3D flow over a swept wing 



b) leading-edge run-back flow 



Figure 10: Schematic diagram for the synthesis of 
the “quasi 3D ” flow over a NACA 0012 swept 
wing. The synthetic flow profile at a given cross 
section (i.e. line aa) mainly contains a 2D BL flow 
over a NACA 0012 airfoil with some prescribed 
spanwise velocity distributions. 
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Figure ll:The most linearly unstable ice roughness 
elements at various chordwise (s*) locations (a) 
s*=0.1mm, (b) s*=4 mm, (c) s*=8 mm, and (d) 
s*=12.6 mm with different film thickness. 
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Figure 12: The growth rates of most linearly unsta- 
ble ice roughness elements at various chordwise 
(s*) locations (a) s* =0.1 mm, (b) s*=4 mm, (c) s*=8 
mm, and (d) s*=12.6 mm with different film thick- 
ness. 
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(1) Conditions from Ref. 17 
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(2) Conditions from Ref. 18 
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Figure 13: Test conditions and results for compar- 
ing the freezing fraction calculation obtained from 
eqn. (18) 17-18 and eqn. (19). 
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